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Abstract 

' We consider real sequences (/„) that satisfy a linear recurrence with 

constant coefficients. We show that the density of the positivity set of 
such a sequence always exists. In the special case where the sequence has 

' no positive dominating characteristic root, we establish that the density is 

positive. Furthermore, we determine the values that can occur as density 
f ^ ' of such a positivity set, both for the special case just mentioned and in 

\ general. 
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1 Introduction and Main Results 

A sequence (f n )n>o of rea l numbers is called a recurrence sequence if it satisfies 
r> ' a linear recurrence 

a- 

fn+h = Clf n+ h-l H h Ch-lfn+1 + Qifn, n>0, 

with constant coefficients Ck € R. (Since we are concerned with questions of pos- 
itivity, we restrict attention to real sequences.) One of the most charming and 
celebrated results in the theory of recurrence sequences is the Skolem-Mahler- 
Lech theorem. It asserts that the zero set 

{n G N : /„ = 0} 
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of a recurrence sequence is the union of a finite set and finitely many arithmetic 
progressions. The recent comprehensive monograph by Everest et al. [3] contains 
references to the substantial literature devoted to this result and to related 
questions. However, not much seems to be known about the positivity set 

{n e N : /„ > 0}. (1) 

In the following section we establish that the density of the set always exists, 
where the (natural) density of a set A C N is defined as 

d(A) := lim x^{n < x : n € A}, 

x — >oo 

provided that the limit exists. 

Theorem 1. Let (/„) be a recurrence sequence. Then the density of the set 
{n G N : /„ > 0} exists. 

Recall that a recurrence sequence can be written as a generalized power sum 

m 

/„ = ]>> fe (n) 7 £, n>0, (2) 

fc=0 

with non-zero polynomials £ (D[n] and roots € C that are roots of the 
characteristic polynomial 

z h - c 1 z h ^ 1 c h _ x z ~ c h 

of the recurrence. The roots of largest modulus are called dominating roots of 
(/„)• 

It does not come as a surprise that recurrence sequences with no positive 
dominating root have oscillating behaviour. Indeed, in section |3 we prove that 
for such a sequence (/„) the densities of Q and the negativity set 

{n e N : /„ < 0} (3) 

are always positive. This generalizes the following known result 00]: If a 
recurrence sequence has at most four dominating roots, and none of them is 
real positive, then the sets and © both have infinitely many elements. 

Theorem 2. Let (/„) be a nonzero recurrence sequence with no positive domi- 
nating characteristic root. Then the sets {n £ N : /„ > 0} and {n € N : /„ < 0} 
have positive density. 

In section 0] we investigate which numbers actually occur as density of the 
positivity set of some recurrence sequence. It turns out that all possible values 
occur, both for sequences with no dominating positive root and in general. 

Finally (section we return to the Skolem-Mahler-Lech theorem. Our 
approach yields the following weak version: The density of the zero set exists 
and is a rational number. 

The conclusion hints at algorithmic aspects of the positivity of recurrence 
sequences. 
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2 The Density of the Positivity Set 



Notation: We write /„ = if /„ = for all n > 0. The Lebesgue measure of a 
set B C R m is denoted by X(B). 

The goal of this section is to prove Theorem ^ Dividing /„ by n |7i| n , 
where 71 is a dominating root of f n and D is the maximal degree of the Pk with 
|7fe| = |7i|, we obtain from JSJl 

d 

n" 13 |7ir n /n = X! ai C0S ( 27T ^ n + ft) + « - r n , 

i=l 

where r n = 0(l/n) is a recurrence sequence, 6\, . . . , 64 are in ]0, 1[ and <Zj, Pi, v G 
R. From now on we will assume w.l.o.g. D = and I71I = 1. As a first step we 
get rid of any integer relations that the 0j's might satisfy. 

Lemma 3. Let 61, ... ,64 be real numbers. Then there is a basis {ri, . . . , r m+ i} 
of the %-module 

M = % + Mi + . . . W d 

such that l/r m+ i is a positive integer and l,7i, ... ,T m are linearly independent 
over (Q. 

Proof. M is finitely generated and torsion free, hence it is free Theorem III. 7.3]. 
Let {ai, . . . , a rn +i\ be a basis. Since 1 £ M, there are integers e%, . . . , e rn +i 
such that 

eiai H h e m+ ia m+ i = 1. 

We complete (ei/g, . . . , e m +i/g), where (7 := gcd(ei, . . . , e m +i), to a unimodular 
integer matrix C with last row (ei/<?, . . . , e m +i/g) [7| §XXI.3]. Then 

(n, . . . , r m+ i) T := C(«i, . . . , a m+ i) T 

yields a basis of M with r m+ i = G (Q. Now suppose 

um + h u m r m = u 

for integers u\, . . . , u m , u. Since u has also the representation 

ugr m+ i = u, 

it follows u\ = ■ ■ ■ = u m = u = 0. □ 

Take ti, . . . , r m+ i as in Lemma |3 with r m+ i — 1/g. Roughly speaking, we 
have put all integer relations among the 0i into the rational basis element T m+ \. 
There are integers 6^ with 

m+l 
Oi = X! b :~r 
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Now we split the sequence (/„) into the subsequences (f gn +k)n>o for < k < g. 
We have 

f gn+k — G n S ni 

where s n := r gn+ k and G n is the dominant part. Defining the integer matrix 

B := {gb ij ) 1 < i < d e TL d * m 

l<j<m 

and the real vector c = (c l7 . . . , q) with 

m+l 

Ci := 2irk °ij T j + Pi, 1 < z < cZ, 

.7=1 

it can be written as (cos is applied component wise) 

G n = a T cos(27rnBr + c) + v. 

We show that the density of {n £ N : f gn +k > 0} exists for each k. Since s n is 
a recurrence sequence with fewer characteristic roots than /„, we may assume 
inductively that d({n £ N : s n < 0}) exists. Thus, if G n is the zero sequence, 
we are done. Now let k be such that G n is not the zero sequence. It is plain 
that G n = H(jit), where 

H (t) := a T cos(27rBt + c) + v, t e [0, l] m . 

The following theorem shows that the function H can be used to evaluate the 
density of the positivity set of G„, which equals, as we will see below, that of 
the set {n : f gn +k > 0}. 

Theorem 4 (Kronecker-Weyl). Let ri,...,r m be real numbers such that 
1, Ti, . . . , r m are linearly independent over Q. Then for every Jordan measurable 
set A C [0, l] m we have 

d({n eN:nT mod 1 e A}) = X(A). 

Proof. We refer to Cassels 2J Theorems IV.I and IV.II]. □ 

We define 

L £ := {n G N : G n > e} and S s := {n € N : |G„| < e}. (4) 
The corresponding sets for the function H are defined as 

L £ := {t S [0, l] m : H(t) > e} and S e := {t e [0, l] m : |if(t)| < e}. 
Since for all e > 

L e = {n S N : nr mod 1 € L e }, 
we have d(L e ) = A(L £ ) for all £ > by Theorem^] Similarly, 

d{S E ) = X(S e ), e > 0. (5) 
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Note that the boundary of the bounded set S e (respectively L e ) is a Lebesgue 
null set (as seen by applying the following lemma with F(t) = H(t) — e), hence 
S £ and L £ are Jordan measurable, and Theorem0]is indeed applicable. Lemma|S] 
seems to be known [H], but we could not find a complete proof in the literature. 

Lemma 5. Let F : M m — > ]R be a real analytic function. Then the zero set of 
F has Lebesgue measure zero, unless F vanishes identically. 

The proof of Lemma [3] is postponed to the end of this section. Since G n is 
not the zero sequence, the function H does not vanish identically on [0, l] m . By 
the Lebesgue dominated convergence theorem and Lemma |5] we thus find 

lim \{S e ) = and lim \{L e ) = \(L ). 

s^O e— »0 

This yields d({n E N : G n > s n }) = \(Lo) by the following lemma, which 
completes the proof of Theorem ^ 

Lemma 6. Let G n and s n be real sequences with s n = o(l) and let L E , S e be 
as in Suppose that d(L £ ) and d(S £ ) exist for all e > 0, and that 

limd(L e ) — d(L ) and lim<i(5 e ) = 0. 

Then 

d{{n e N : G n > s n }) = d(L ). 
Proof. For any set A C N we write A(x) :~ {n < x : n <E A}. Define 

P := {n e N : G n > s n }. 
Let e > be arbitrary. Take no such that \s n \ < e for n > uq. It follows 



tt-P(x) = (t{n < n : G n > s n } + jj{n < n < x : G n > e} 
+ U{n < n < x : s n < G n < s}, 

hence 

\iP(x)-tL e {x)\<$S e (x)+o(x) 
as x — > oo. Thus we have 

\x-HP(x) ~ d(Lo)\ < \x-HP{x) ~ x~HL s {x)\ + \x- x $L e (x) - d(L )\ 
< x-^S e (x) + \x-^L £ (x) - d{L )\ + o(l). 

The right hand side goes to 

d(S s ) + \d(L e )-d(L )\ 

as x — > oo. By assumption, this can be made arbitrarily small, which implies 
d(P)=d(L ). ' □ 
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Proof of Lemma\^ For m = 1 this is clear, since then the zero set is countable. 
Now assume that we have established the result for 1, . . . , m — 1. Put 

V := {(t 2 , t m ) E R 1 ™" 1 : F(-, t 2 , . . ■ , t m ) vanishes identically }. 

Take a real number s such that F(s, ■,...,■) is not identically zero. Clearly, 
F(s, *2, ■ • ■ , t m ) = for all (i 2 , ■ • ■ , t m ) S V. By the induction hypothesis, this 
implies X(V) — 0. Note that V is closed, hence measurable. Since F is real 
analytic in the first argument, we have 

/ Xz(ti,...,t m )d\(h) = 

JR 

for all fa, ■ ■ ■ , t m ) $ V, where \z is the characteristic function of the zero set 

Z := {(ti, ,.,i m )6l m ; F(ti, . . . ,t m ) = 0}. 
Since V has measure zero, this implies 

f ... f xz(h,...,t m )d\(t 1 )...dX(t m ) = o. 

This argument works for any order of integration, hence we obtain J" Rm xz = 
by Tonelli's theorem. □ 

3 Sequences with no Positive Dominating Root 

In this section we prove Theorem^ We begin by settling the special cases where 
the 6i are all irrational or all rational, and then put them together. 

Lemma 7. Let 9i, . . . ,9<i be irrational numbers, and let ai, Pi be real numbers 
such that the sequence 

d 

u n = ai cos(27r#i7i + /3j) 

is not identically zero. Let further r n be a recurrence sequence with r n = o(l). 
Then the set {n € N : u n > r n } has positive density. 

Proof. Proceeding as in the proof of Theorem ^ we can write 

G n := Ug n+ k = a T cos(27rnBx + c), 

where B is an integer matrix no row of which is zero, c is a real vector and 
1, Ti, . . . , t t „ are linearly independent over (Q. If k is such that G n = u gn +k = 0, 
then the density of {n € N : G n > s n }, where s n — r qn+ k, exists by Theorem^ 
but may be zero. Now choose a fco such that the corresponding sequence G n — 
Ugn-\-k is n °t the zero sequence. We have G n — H(jit), where 

H(t) := a T cos(27rBt + c). 
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Moreover, with the notation of the proof of Theorem we have 

d({n € N : G n > s n }) = A(Z ). 
The function H is not identically zero on [0, l] m . But 



f ... f H{t 1 ,...,t m )dt 1 ---dt m = 0, 
Jo Jo 



(6) 



because no row of B is the zero vector. Hence H has a positive value on [0, l] m , 
and since it is continuous, we have X(Lq) > 0. □ 

Observe that the integral in © need not vanish if B has a zero row, which 
can only happen if the 9i corresponding to this row is a rational number. This 
is the reason why we consider rational 9i's separately. 

Lemma 8. Let 9i,...,9d be rational numbers in ]0, 1[, and let ai,Pi be real 
numbers such that the purely periodic sequence 



= 22 a i cos(27r#iro + ft) 



i=i 

is not identically zero. Then u n has a positive and a negative value. 
Proof. By the identity 

cos 2Z^E + i sin = ^ ^ ikp/q = 0, 

fc=0 k=0 k=Q 

valid for integers < p < q, and the addition formula of cos we obtain 

u -\ 1- U q -l = 0, 

where q is a common denominator of 9% , . . . , 6^. □ 
Proof of Theorem^ It suffices to consider the positivity set. We may write 

fn ~l~ 

where r n = o(l) is a recurrence sequence, 

d 

u n = ^^ai cos(27r^n + ft), 

i=l 

e 

v n = ^ ai cos (2ir9in + ft), 

i=d+l 

9\ 1 . . . ,9d are irrational, 9d+i, ■ ■ ■ ,9 e are rational numbers in ]0, 1[ with common 
denominator q > and u n + v n ^ 0. If v n = 0, then the result follows from 
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Lemma [7] Now suppose v n ^ 0. Then for each k the density of the set {n £ 
IN : fqn+k > 0} exists by Theorem ^ By Lemma |H1 there is ko such that 
v qn+k„ — v > 0. It suffices to show that the set {n G N : f qn +k a > 0} has 
positive density. This is clear if u qn+ko = 0. Otherwise, notice that 

{n € N : f q n+k > 0} 3 {n G N : u g „+fco > r qn+ko }, 
and the latter set has positive density by Lemma [7| □ 



4 The Possible Values of the Density 

In this section we investigate which values from [0, 1] occur as density of some 
recurrence sequence. In its basic form, the question is readily answered: 

Example 9. Let w be a real number and define 

f n := sin(27rnV / 2) — w. 

Then, by Theorem^ 

d({n € N : /„ > 0}) = \{{t S [0, 1] : sin(2vrf) > w}) 

1 w < -1 

i — — arcsin w — 1 < w < 1 . 

w > 1 

Since the range of arcsin is [— -f , /or every k € [0, 1] i/i«s yields a recurrence 
sequence f n such that 

d{{n g N : /„ > 0}) = k. 

The following proposition generalizes this example. Note that the density 
of the zero set of a recurrence sequence is always a rational number by the 
Skolem-Mahler-Lech theorem. 

Proposition 10. Let n be a real number and r be a rational number with 
< n,r < 1 and K + r < 1. Then there is a recurrence sequence (/„) such that 

d({n € N : /» > 0}) = re and S N : /„ = 0}) = r. 

Proof. Suppose that r — p/q for positive integers p and g. As seen in Example^! 
there is a recurrence sequence (g n ) such that the density of the zero set of (g n ) 
is zero and the density of its positivity set is — r) (The case r = 1 is trivial). 
The interlacing sequence 



fbn+k 



< k <p 
g n p< k < q 



8 



is a recurrence sequence |31 section 4.1]. Clearly, the density of its zero set is r, 
and the density of its positivity set is 

d({n e N : /„ > 0}) = x JL. = K , 

q 1 — r 

as required. □ 

If we restrict attention to sequences without dominating real positive roots, 
then Theorem |3 tells us that the density of the positivity set can be neither zero 
nor one. Still, all values in between occur. 

Theorem 11. Let n G]0, 1[. Then there is a recurrence sequence (/„) with no 
positive dominating characteristic root and d({n € N : /„ > 0}) = k. 

Proof. Let e > be arbitrary. We define a function H on [0, |] by 

f (e-1) 2 



H{t) 



(i-f) 0<t<§ 
e-2t §<*<§ 



and extend it to an even, 1-periodic function H on R (see Figure It is 
continuous and satisfies 

j H(t)dt = and \({t e [0, 1] : H{t) > 0}) = e. 
Jo 



1 




—+6-2 




e 










-1 +e 





Figure 1: The function H 

Expanding H into a Fourier series, we find that there are real a,j such that 
H is the pointwise limit of 

m 

H m {t) := ^2 a o cos(27rji) 

3=1 
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as m — ► oo. Since the zero set of H is a null set, the Lebesgue dominated 
convergence theorem yields 

lim A({ie [0,1] : H m (t) >0}) = e. 

m — >oc 

Wc fix an m such that 

\{{t£ [0,1] : H m (t) > 0}) < 2e. 

The function 

cj ) (A 1 ,... 1 A m ) :=A({tG [0,l]:]TA,-coB(27rjt)>0}) 

i=i 

is continuous on R m \{0}. To see this, observe that is continuous at all 
points (Ai,...,A m ) for which 2j=i -4j cos(2wjt) is not identically zero and 
appeal to the uniqueness of the Fourier expansion. Since 0(1, 0, . . . , 0) = h 
and 4>{a\ 1 . . . , a m ) < 2e, the function <^> assumes every value from [2e, i] by the 
intermediate value theorem. 

Hence the positivity sets of the sequences 

m 

fn ■= Y2 A 3 c °s(2TrjnV2) 

assume all densities from [2e, \] for appropriate choices of (Ay, . . . , A m ) by The- 
orem0] Repeating the whole argument with —H instead of H yields the desired 
result for k G [|, 1 — 2s]. Since e was arbitrary, the theorem is proved. □ 

5 A Weak Version of Skolem-Mahler-Lech 

Without using the Skolem-Mahler-Lech theorem, it follows from Theorem ^ 
that the density of the zero set of a recurrence sequence (/„) exists. We can 
show a bit more with our approach. Recall, however, that we only deal with 
real sequences, whereas the Skolem-Mahler-Lech theorem holds for any field of 
characteristic zero. 

Proposition 12. The density of the zero set of a (real) recurrence sequence 
(f n ) is a rational number. 

Proof. Let k be a natural number, and let g, G n and s n be as in the proof of 
Theorem ^ If k is such that G n = 0, then the density of the zero set of f gn +k 
is rational, since we may assume inductively that the density of {n : s n = 0} is 
rational. 

Now suppose G n ^ 0. The zero set of f gn +k can be partitioned as 

{n e N : G n = s n } = {n : G n = s n , |G„| < e} U {n : G n = s n , \G n \ > e}, 
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where e > is arbitrary. The latter set is finite, and the first one is contained 
in S £ , defined in (@J. Hence 

d({n g N : G n = s n }) < d{S E ) 

for all £ > 0. But we know that lim e ^o d(S £ ) — from the proof of Theorem^ 
which yields 

d({n £ N : G„ = s„}) = 0. 

Thus, the zero sets of all subsequences (f gn +k)n>0: < k < g, have rational 
density, which proves the desired result. □ 

6 Conclusion 

There is no algorithm known for deciding whether /„ > for all n, nor has 
the problem been shown to be undecidable. When we are talking about algo- 
rithmics, it is natural to assume that the recurrence coefficients and the initial 
values are rational numbers. In this case Gourdon and Salvy [H] have proposed 
an efficient method for ordering the characteristic roots w.r.t. to their modu- 
lus. Thus, the dominating characteristic roots can be identified algorithmically. 
If none of them is real positive, then we know that the sequence oscillates by 
Theorem On the other hand, sequences where a positive dominating root is 
accompanied by complex dominating roots seem to pose difficult Diophantine 
problems. For instance, we do not know if the sequence 

/„ :=cos(27T0n) + l + (-i) n (7) 

is positive for 9 = v2 ; say. It can be shown, however, that the set of #'s for 
which the corresponding sequence (/„) (defined by 10) is positive has measure 
zero Theorem 7.2]. 
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